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Two dimensional suspensions of spherical colloids subject to periodic external fields exhibit a rich 
variety of molecular crystalline phases. We study in simulations the ground state configurations 
of dimeric and trimeric systems, that are realized on square and triangular lattices, when either 
two or three macroions are trapped in each external potential minimum. Bipartite orders of the 
checkerboard or stripe types are reported together with more complex quadripartite orderings, and 
the shortcomings of envisioning the colloids gathered in a single potential minimum as a composite 
i-O i rigid object are discussed. This work also sheds light on simplifying assumptions underlying previous 

^ ■ theoretical treatments and that made possible the mapping onto spin models. 
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I. INTRODUCTION 

o : 

Colloidal suspensions provide valuable systems for the study of collective effects and phase transitions (see e.g [l[ or 
more specifically Q and references therein). Recent advances in optical trapping techniques have enlarged the whole 
field and opened the possibility to obtain quasi two-dimensional systems, which can furthermore be subject to an 
I , external potential. A possibility to realize such external perturbations is through the interference of laser beams, the 
typically micron-sized colloids being attracted to the regions of highest intensity: for instance.unusual phases have 
been reported in the case of ID troughs [H, such as floating solids or locked smectic phases [3, HI- 

The situation of a two-dimensional periodic substrate, some instances of which can be realized by 3 laser beams, 
is remarkably rich and has been addressed experimentally @ , numerically @, H, Q and analytically [1(1 HH ■ The 
control parameters governing the static behaviour are numerous : filling fraction (mean number of colloids per 
substrate minimum), pinning amplitude (trap strength, increasing with laser intensity), temperature, and concerning 
the substrate geometry, lattice spacing and aspect ratio (i.e rectangular unit cell instead of a square one). In addition, 
the colloids considered are highly charged objects, so that the Debye length of the suspension, modified by changing 
the salt content, is a crucial parameter. In the present paper, we concentrate on the experimentally relevant ground 
state of those 2D systems, where the long range orientational order observed has been coined "colloidal molecular 
crystal" . If the external laser potential is strong enough, the colloids are irreversibly bound to the potential minima. 
In cases when more than one colloid is trapped in a single minimum we speak of a "colloidal molecule". Its size is 
00 ' determined by the interplay between light forces and interparticle repulsion. Figure [TJa) provides an illustration of 
typical phases observed when the ratio of the total number of colloids in the system to the total number of external 
potential minima is exactly two. Whereas defects are present at finite temperature, we observe that the ground state 
is free of such objects, so that each trap captures exactly two colloids. 

The purpose of the present paper is twofold. The first goal is to provide a thorough numerical investigation 
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of orientational ordering in colloidal molecular crystals on two-dimensional periodic substrates. We will focus on 
stoichiometries 2 and 3, where dimers and trimers are formed, respectively, with the underlying square or triangular 
symmetry for the potential. We will also address the situations where the corresponding lattice unit cell is distorted 
by changing the aspect ratio a. The second objective is to critically test several assumptions that led to the theoretical 
frameworks used in Refs. [HI [HI to study such problems. In section UH we will start with an approach, common 
to Refs fill [Tl| . where the composite objects (n-mers) formed in each trap are considered as rigid entities with an 
orientational degree of freedom only, while the confining potential is taken into account implicitly. This allows for 
a significant reduction of the complexity of the problem. Reference phase behaviour is thus obtained, which will 
be tested against more realistic simulations in Section [TTTT where each colloid is resolved and the confining potential 
explicitly accounted for. Finally, we will also test the relevance of the various effective potentials used in to 



'Electronic address: samir.elshawish@ijs.si 



2 



construct a tractable Hamiltonian allowing for analytical progress. Our main finding and the ensuing consequences 
will be summarised in section |Vl 




FIG. 1: a) Several representative orientational structures observed on a 2D rectangular lattice, for a stoichiometry n — 2 
(dimers). For visualisation purposes, the two colloids forming a dimer are hereafter linked by a thick line. The parameters 
associated with the different phases are the following: for the structure denoted Po,o (ferromagnetic): nd = 0.25, kI — 6, a — 0.9; 
Pe,e (tilted ferromagnetic): nd = 1, nl = 4.5, a = 1; Po 1 ,e 2 : nd — 1.5, kI = 6, a = 1; and for the Po,tt/2 (antiferromagnetic) 
phase: nd — 0.5, kI = 6, a = 1. 1/k is the inverse Debye length and represents the range of the electrostatic screened 
interactions between colloids, b) Schematic drawing showing dimers in the traps (shaded circles). I is the distance between the 
neighbouring traps and 2d is the size of the "colloidal molecule". 



II. THE RIGID n-MER APPROACH 

The integer number n of colloids that gather in a light potential minimum are subject to gradient forces arising 
from the dielectric mismatch between colloids and the solvent, and light pressure [l3j . Under such forces alone, the 
colloids have a preference for regions of highest laser intensity. In addition, the (spherical) colloids interact through 
strong mutual repulsion, considered to be of a screened Coulomb form [l4| 

$c = ifE exp(-^) 

where is the distance between the centres of masses of macroions i and j, 1/k is the Debye length [l5l |. and K is an 
irrelevant prefactor. To a large extent, colloids in a given trap experience the repulsion from their n — 1 "trap-mates" 
only, and interactions with colloids in other traps are of little relevance to determine the shape that the n-mer adopts 
within a trap. We therefore assume here that the antagonistic effects of light interaction and Coulomb repulsion lead 
to the formation of rigid composite objects with n-fold symmetry (e.g. equilateral triangles for trimers when there are 
3 colloids per trap). This is the basic assumption of Refs [lfj. for which the relevant ground state dimensionless 
parameters are kI, nd and aspect ratio a, where I is the distance between the centres of two adjacent traps, d is the 
distance between the centre of the trap and the centre of one of the colloids forming an n-mer (so that in the dimeric 
case, 2d is the dimer extension, see Fig. [Ub)). 

The rigidity assumption lumps substrate potential effects into the length scale d (increasing the well's amplitude 
decreases d). The position of each n-mer is then described by a unique angular coordinate, while its centre-of-mass 
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always coincides with trap minimum. Unlike in Ref. [ll| , we do not assume that such angles are discrete with values 
dictated by the lattice geometry: it is indeed of interest to realize that Coulomb repulsion alone is able to select well 
defined orientations that do not match any of the underlying lattice principal direction. Finding the ground state of 
a system of iV n-mers therefore amounts to minimising the energy function (TTJ) with respect to N angles fl6| . 

For the simplest case of dimers (n = 2) on a square (a — 1) or rectangular (a ^ 1) lattice, we found that in the 
lowest energy configuration -obtained by a simulated annealing technique , the system adopts a bi-partite structure 
of the checkerboard type. The long range orientational order is therefore characterised by two angles 8\ and 82, which 
define a phase denoted Pg lt g 2 , see Fig. [TJ Interestingly, the angles 8\ and 82 are constant in some parameter range, 
while they vary continuously in other parameter regions. This is illustrated in Fig. [2j The low nd regions of graphs a) 
and b) are such that (81,82) = (0,7r/2), a situation coined "antiferromagnetic" in previous studies and illustrated in 
the lower right corner of Fig. [TJ Starting from nl = 6 with a square lattice (Fig. Ob) and slightly distorting the lattice 
into a rectangular one with aspect ratio 0.9, we observe that the antiferromagnetic Pq^/2 phase disappears and turns 
into a ferromagnetic one Po,o, for small enough dimer extensions (small Kd), see Figs. [2]-d) and -e). For higher Kd, a 
tilted ferromagnetic Pg,e phase appears the most stable, and changes into a tilted antiferromagnetic Pe 1 ,e 2 (Figsl2]-d, 
-e). Mild parameter differences therefore trigger significant orientational changes: for instance, at kI — 6, all phases 
are antiferromagnetic-like (see graph b)) while at kI = 4.5 (graph a)) there exists a window around kg? = 1 with a 
ferromagnetic ordering. Note also that the upper part of phase diagram c) corresponds to cases where d > 1/2, which 
is ruled out since a given trap cannot extend further than half the inter-trap distance. Similarly, the forbidden upper 
region of diagram f) corresponds to d > al/2. 
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FIG. 2: Ground state of dimers (n = 2) on a square or rectangular lattice, as obtained from simulated annealing. Graphs a) 
b) d) and e) exhibit the dependence of characteristic angles 0i and 62 on dimer spacing d, lattice constant I and aspect ratio 
a. The first row is for a square lattice (a = 1) with (a) kI — 4.5 and (b) kI — 6. The corresponding phase diagram is shown in 
(c). The second row is for nl — 6 and (d) a = 0.9, (e) a = 0.7, while (f) shows the phase diagram. As in subsequent figures, 
the symbols refer to parameters used for the more complete Monte Carlo simulations of section IIHI which lift the rigidity 
assumption. Circles: "rigid" simulated annealing and "flexible" Monte Carlo simulations agree; triangles for parameters where 
they do not agree. Here, no triangles are reported since both approaches yield similar results, see below. 

When the lattice geometry is triangular, we have observed the formation of stripes @, 0], and the corresponding 
phases are denoted Sg lt g 2 (the stripes are of the same bi-partite family as the order reported on the square lattice, 
since the checkerboard structure itself is made up of parallel stripes). The main results are summarised in Fig. [3] 
The upper left and lower left regions of graphs a) and b) respectively are forbidden regions where some configurations 
of the rigid trimers in neighbouring traps would lead to overlaps. For a given aspect ratio, two ground states are 
generically observed : a herringbone order (see e.g. the upper inset of Fig. [3]-a)) and a ferromagnetic one where dimers 
align. The sequence of these two phases as parameters are modified is quite complex, particularly so when the aspect 



4 



ratio a is changed. The lattice with maximal symmetry a = 1 appears singular in that the only phase selected there 
is the herringbone one. This qualitatively confirms the results of The case of trimers on a square or rectangular 
lattice is somewhat simpler, with only the stripe order observed (see Fig. ffl-a)). When the trimers are put on 

the triangular lattice, several stripe phases are possible, but roughly speaking, there is a unique type of order for a 
given aspect ratio, see Fig. HJ-b). 




Kd a 

FIG. 3: Phase diagram for rigid dimers on a triangular lattice for (a) a — 0.9 and (b) nd = 1. The insets show the representative 
stripe configurations. The parameter dependence of bipartite angles along the td = 6 line of phase diagrams (a) and (b) is 
shown in (c) and (d), respectively. Circles have the same meaning as in Fig[2] 




FIG. 4: Phase diagram of rigid trimers on (a) rectangular lattice for a — 1 and (b) triangular lattice for nd — 1. The stripe 
phases are sketched in the insets, where the centre of masses of the three colloids gathered in a potential minimum are linked 
by a thick line. Here again, circles indicate an agreement between "rigid" and "flexible" simulations, and triangles report a 
mismatch. 



III. TOWARDS A MORE COMPLETE DESCRIPTION : EFFECTS OF FLEXIBILITY 



To test the relevance of the rigidity assumption we now explicitly take into account the confining potential in the 
simulations. In order not to introduce any orientational bias, we consider an isotropic (harmonic) confining potential 
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®L(r) — Vo(Kr) 2 , so that the total dimensionless energy of colloid i reads: 

ei = Ek/K = A(KSi) 2 + V CXP( ~^ j) , (2) 

. . . K>Tij 

where Si denotes the distance between the centre of mass of colloid i and the trap centre. The relevant parameters 
are now hd, A = Vo/K which measures the relative strength of the light confinement against the Coulomb repulsion 
and therefore sets the n-mer size denoted d in section [II] and again the aspect ratio. We have performed Monte Carlo 
simulations of this "flexible" model, which has 2nN degrees of freedom whereas the "rigid model" only has N degrees 
of freedom. 

The results are shown in Figs. 0(5] and [71 and depend on the particular situation. The comparison flexible/rigid 
was performed by first implementing simulated annealing for the flexible model, measuring the resulting n-mer size, 
and using the corresponding value of d in a rigid-model simulation. For dimers on the square/rectangular lattice, we 
always found an excellent agreement between both routes, see Fig. [5] where three typical bipartite orders are shown 
in the insets. On the triangular lattice, we found that the qualitative features put forward in section [IT] remain correct 
(see e.g. the configuration shown in the inset of Figure , while on closer inspection, some differences arise. Indeed, 
the histograms of distances to trap centre S and tilt angle 8 in Figure [HI clearly reveals the mismatch between the 
rigid and flexible ground states that are respectively of bipartite and of quadripartite type (the inset shows the new 
unit cell by a shaded area). Correlating the histograms to the snapshot of the inset, it appears that the top row of 
the inset corresponds to the most intense peak in the P{9) distribution, and to the two extreme peaks in P(S) at 
6 ~ 0.65c? and 1.35(1 Conversely, the second row from the top has colloids contributing to the two smaller peaks in 
P(0) that weight half the previous one, and contributing to the peaks at 8 ~ d. It can be seen that the average 6 is 
d, as it should from the definition of d in the flexible case. 
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FIG. 5: Comparison between the ground states for dimers on a square substrate, obtained considering the full "flexible" model 
including explicitly the trapping potential, against the restricted approach where the n-mers are considered as rigid objects. 
The curve corresponds to the rigid scenario, while the circles are for its flexible counterpart. Here, nl — 4.5 and (a) nd = 0.4, 
(b) Kd = 0.97, and (c) Kd = 1.62. 

Figure [JJ shows similar results for trimers on both rectangular and triangular lattices. On the rectangular lattice, 
the main difference with the rigid case (where the trimer centre is imposed to coincide with the trap centre) , lies in an 
off-centre shift, see the inset of graph a). Graphs a) and b) corroborate and quantify this visual observation. On the 
triangular lattice, the effect is more spectacular. While the rigid model leads to an order where all trimers align in 
the same direction (<So,o fashion), due account of the internal flexibility of the trimers yields an order that is visually 
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FIG. 6: Emergence of tetrapartite ground-state configurations for flexible dimers on a triangular lattice with kI = 6, nd = 1.94, 
a — 1. The "rigid" configuration is shown for comparison with dimers indicated by lines, while the discs materialise the colloids 
in the flexible ground state. Graph a) shows the probability distribution function of colloid position in the trap (5 = for the 
trap centre). Graph b) is for the distribution function of the angle between a reference direction and the line joining the trap 
centre to a given colloid. As for the inset showing the configuration, the flexible data are shown with discs and the rigid ones 
are shown by the continuous line. 



reminiscent of a stripe phase , but that is in reality more complex. This can be appreciated by comparing the 

configurations in the first and third rows from the top of the inset in graph c) (or equivalently second and fourth): 
an off-centre shift of the triangles may be observed, and substantiated by the histograms of graphs c) and d). Here 
again, the order selected is of tetrapartite form, with a unit cell shown by the shaded area, and not bipartite as found 
in section HU 
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FIG. 7: Same as Figure[6]for trimers on a rectangular lattice (graphs a) and b) for which nd = 1.88, nl = 6 and a = 0.9), and 
for trimers on a triangular lattice (graphs c) and d) for which nd = 1.66, kI — 6 and a = 0.75). 
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IV. ON THE RELEVANCE OF DIFFERENT APPROXIMATION SCHEMES 

In the previous section, we have shown that assuming the trapped colloids to form a rigid composite object may be 
incorrect in some cases. Here, we address the applicability of two other classes of simplification of the original model 
provided by the energy function @. 



A. The nearest neighbour interactions 



The first simplification, considered in both Refs. [10( and consists in restricting colloid interactions to partners 
in nearest neighbour traps. Given the exponential character of screened Coulombic law ([I]), this seems an a priori 
reasonable assumption, provided both the distance between adjacent traps (I) and the closest distance between colloids 
(/ — 2d) are large enough compared to Debye length 1/k. The comparison between the two approaches is displayed in 
Fig. [5] for dimers on a square or rectangular lattice, and on Fig. [5] for dimers on a triangular lattice. We consider only 
the rigid dimers, which has been shown to be sufficient in section HTI1 see Fig. [5]). The above argument leads to believe 
that for large enough kI and small enough nd, interactions beyond the nearest neighbours should be immaterial. This 
can be observed in Fig. [8l where the angles &\ and 02 (those of the bipartite checkerboard structure) are plotted 
against nd for six different (kI,ci) combinations. The nearest neighbours approximation generally works well when 
nd <C nl with, however, the surprise that at nl — 4.5 and small nd, the nearest neighbour route is quantitatively wrong, 
leading to a tilted ferromagnetic phase P^u^/a instead the antiferromagnetic Po, n /2 ( see Fig- 130)- A comparison 
between the energies of the different phases is performed in Fig. 1101 Surprisingly, the relative energy difference 
between ferromagnetic and antiferromagnetic phases is minute (see the y-scale) and explains why including second 
nearest neighbours is essential to get the correct phase behaviour. The dashed line in Fig. [TO] is obtained by summing 
over the four nearest neighbour traps and over the four next-nearest corner traps in the rectangular array. Restricting 
to nearest neighbours leads to the prevalence of the tilted ferromagnetic (symbols), which is incorrect, although 
energetically very close to the antiferromagnet. When kI is increased, this quasi degeneracy disappears, and including 
next nearest neighbours in the analysis becomes irrelevant. 
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FIG. 8: Comparison of the nearest neighbour assumption (shown with symbols) and the correct result (shown with solid lines) 
for dimers on a rectangular lattice. The top row is for a square geometry (a — 1), with a) kI = 4.5, b) kI = 6 and c) nl = 8. 
Second row is for nl = 6, tuning the aspect ratio: d) a — 0.9, e) a = 0.85 and f) a = 0.7. 



In the case of the triangular lattice, the nearest neighbour approximation is more efficient than on the rectangular 
lattice, as can be seen in Fig. [SJ The reason for this is most likely that in the triangular lattice where there are six 
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FIG. 9: Comparison of the nearest neighbour assumption (shown with symbols) and the correct result (shown with solid lines) 
for dimers on a triangular lattice. The parameters are the same as in Fig. [S] 




FIG. 10: Comparison of the energies of the phase (denoted AF) and the phase (denoted F*) for dimers on a 

square lattice (a — 1) with kI = 4.5 and small ltd. This parameter range corresponds to the small nd region of Fig. [8]-a). The 
symbols are for nearest neighbour interactions only, the dashed line adds the next shell of nearest neighbours (the four corner 
sites in the square), while the solid curve is for the case where all neighbours in all traps are considered. Due to the extremely 
small energy differences Vaf — Vf* summation over more neighbour shells is required and the nearest neighbour approximation, 
although intuitively expected to be valid, fails. 



nearest neighbour sites, the next-nearest neighbour distance is y/3 times larger than the nearest neighbour one. On a 
square lattice on the other hand, there are four nearest neighbour sites and the next-nearest distance is only a factor 
V2 larger. The situation is thus less favourable for nearest neighbour truncation. 
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B. The large-distance approximation 



The second simplification enforced in [Io| with the purpose to allow for the construction of a tractable Hamiltonian, 
which can be mapped onto a spin model, amounts -within the rigid scenario- to considering the leading order term 
only in the large distance expansion of the interaction potential. The reason for doing so is again to focus on the 
large td case, with the technical bonus that then, interactions between 2 n-mers in different traps may be written in 
a factorized way. To be specific, the large distance potential of interaction between two dimers labelled 1 and 2, of 
size 2d, with centre-of-mass/centre-of-mass separation r reads 



V\2 = cosh[ft<icos(#i2)] cosh \nd cos (#21) 



(3) 



where r = |r| and 6ij is the angle between vector r and the direction defined by dimer i. Similar considerations prove 
generically fruitful to discuss interactions between anisotropic colloids in the low density regime [13, EH . It can be seen 
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FIG. 11: Effect of considering the large distance approximation (LDA) of the interaction potential between dimers (on rectan- 
gular lattices) on the bipartite angular ordering. The parameters are the same as in Fig. [8] The truncated potential results 
with expression ([3]) are shown with the symbols: triangles restricting to nearest neighbours, and circles including all neighbours. 
The results are compared to the "all neighbours" simulation of the rigid dimers (solid lines). 



in Figure QT] that for low nd, the corresponding predictions fare favourably against the results of the full expression 
(JTJ. However, the truncated approach always predicts the antiferromagnetic phase in the square case 0, and misses 
the tilted antiferromagnet (see the top row). On a rectangular lattice, it correctly captures the transition from a 
ferromagnet Po.o to tilted ferromagnet, but overestimates the threshold (see the bottom row), to such an extent that 
it can exceed al/2, the maximum allowed value of d (see graph f)). The failure at larger nd has a different origin from 
the one observed in Fig. [8j and could have been anticipated by computing the sub-dominant terms in expression (|3j). 
For the sake of clarity, we focus here on the electric potential created by a dimer, a simpler but related object than 
the dimer-dimer potential considered in {SJ; the sum of the two screened Coulomb terms associated to each colloid 
may be written 

find) 2 \ e~ Kr 

V(r) ~ cosh[/cdcos(0)l exp V^(cos 2 0-1) , (4) 

\ 2nr J nr 

where 9, defined as above, is the angle under which the dimer is seen from a distance r. Neglecting the second factor 
on the right hand side (which leads to Eq. © for the interaction) is only justified provided (nd) 2 < td. Finally, we 
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note that although the parameters in Figs. [8l and [TT1 are exactly the same, truncating colloid interactions to nearest 
neighbour traps only makes little difference if the potential is of the form ([3]). This is at variance with what can be 
observed in Fig. [H We are back here to the message conveyed by Figure [lOl that minor modifications of the original 
problem may significantly alter the preferred phase. 



C. The discrete angles approximation 

In reference ill], an effective Potts-like Hamiltonian was constructed by assuming rigid molecules of fixed sizes 
placed on the lattice points, restricting the orientations 6i of the molecules to discrete values compatible with the 
lattice symmetry and by considering the nearest neighbour interactions - Yukawa interactions between colloids in 
neighbouring traps. Since the phase behaviour studied in [ill ] pertains to triangular substrate potentials, the nearest 
neighbour approximation seems well justified, see Fig. [9J The restriction to discrete angles, which rests on the remark 
that the confinement potential exhibits some preferred directions, however deserves a more careful discussion. 

For concreteness, we consider in the triangular geometry case, a confinement potential of the same form as in Ref. 

0: 



= const 



2tt + cos 47r — = + cos 2w 



(5) 



The iso-0L lines are isotropic (circular) in the vicinity of the minima, with an anisotropy that increases with 
increasing energy <f>L. As a consequence, the confinement potential anisotropy is all the more important as the n-mer 
is more extended (large d), but this parameter range corresponds to a situation of strong Coulombic repulsion where 
details of the light potential may not matter. It is therefore not straightforward to anticipate the relative ranges of 
applicability of the discrete angle approximation of [Tlj , and of the isotropic confinement potential approach followed 
in this paper, which leads to Figs. fTlfTTl 

To answer this question, we have minimised the total energy e = A'Y^ li <\>L + <i>Ci with respect to the angular 
orientation of the dimer. Such an approach takes due account of the interplay between the realistic potential ([5]) and 
Coulombic interactions. The ground state results are shown in Fig. Q21 for the same parameters as in the experiments 
of Ref [||: a = 1, nl = 20. The dimers are assumed to form the herringbone structure Sg li g 2 , where due to the 
symmetry, we may assume 02 = n — 9i (see Fig. [3]). For a given ratio A between the light confinement and the 
strength of the Yukawa interaction, the total energy e has been minimised with respect to the radial (red) and angular 
(6?i ) positions of the colloid inside the trap. The resulting angles 6q are plotted against red in Fig fl2T a). see the solid 
line. For comparison, the lower dot-dashed curve is the result obtained with the isotropic potential (angles selected 
by the Yukawa interactions). The upper dashed curve is the constant angle Q\ — 7r/3 assumed in fill ]. It appears that 
the angles are compatible with the pure Yukawa prediction at large values of red. In essence, although the confinement 
potential is more anisotropic at large d (corresponding to small A) , the Yukawa term nevertheless dominates the total 
energy so that the isotropic potential approach becomes correct. Conversely, the discrete angle approach appears 
more relevant for small colloidal molecules. 

In Fig. [T2Tb) we show the radial dependence of the total energy, fixing 8 to its optimal value. The total energy 
E = Ke is obtained assuming the experimentally realistic value K = 10 5 fcsT, the value used also in [ill ]. Larger 
molecules have weaker radial binding energies and are therefore more prone to fluctuate at finite temperature. In 
FigfT2Tc1. we plot the angular dependence of the total energy at fixed red. The depth of the minimum can be 
understood as the angular binding energy. This binding energy is of several ksT at large red, while it is very small 
at small red, in the 9 = n/3 regime [19j . The assumption of discrete angles therefore fails at large red (weak light 
potentials) and leads to the correct results only at small red, where, however, due to the shallow binding in the angular 
direction, there is no a-priori reason to justify it but for ground state properties. 

The parameter regime of the experiments Q is marked as a shaded area in Fig. [T27 a) . Interestingly, the experimental 
situation is somewhat intermediate between both limiting regimes and neither of the approaches seems to be perfectly 
appropriate. However, other realizations of the light confinement are possible, like creating isotropic optical point 
traps instead of using the interference patterns as in [6j. For such an experiment, the present approach - the lower 
curve in Fig |12f a). corresponding to the isotropic trap - would be relevant. 

Finally, we also note that numerical simulations of colloidal molecules in the anisotropic light confinement (J5]) 
have been performed in The parameters studied there (re/ = 2) are nevertheless quite distant from those of the 
experiments Q where nl rs 20. In an aqueous solution, nl = 2 cannot be achieved with 3 nm size colloids as used 
in [||, since the inverse Debye length re is bounded from below by the solvent dissociation. Moreover, in the regime 
nl w 2, the pair-wise additivity assumption leading to the Yukawa potential is questionable, and many-body effects 
could play an important role. 
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FIG. 12: a) The optimal orientation 8 of a dimer in the triangular light potential obtained through minimization of the total 
energy (confinement+Coulombic). The dot-dashed line is for an isotropic parabolic light potential where only the Yukawa part 
of the energy determines the angle 6. The upper dashed line is the constant value 9 = n/3 compatible with the lattice symmetry 
of the light potential. The solid line is for the case of anisotropic confining potential ([5]) where orientations favoured by Yukawa 
interactions and those favoured by the lattice symmetry, compete. The shaded area represents the parameter regime of the 
experiments [f|. b) Radial cut through the total energy landscape for different points in graph a) corresponding to different 
values of the ratio A. The depth of the minimum represents the binding strength of the dimer in radial direction. The sizes 
of the molecule are determined from the location of this minimum, c) Angular cut through the total energy landscape for 
the same points as in graph b). The depth of a minimum represents the angular binding strength of the dimer. This angular 
binding energy is of the order of a few ksT at large nd (point D) and very small for smaller molecules (curves A, B, C). The 
total energy E = Ke in b) and c) is obtained assuming K = 10 5 ksT. 



V. CONCLUSION 



We have provided a numerical analysis of the type of long range orientational orders selected in so-called colloidal 
molecular crystals, where a given integer number n of colloids is trapped in each potential minimum of a light 
lattice. We investigated the cases of dimers (n — 2) and trimers (n = 3) on both square and triangular lattices, 
together with deformed geometries obtained when a given direction of the original lattice is expanded or shrunk 
by a factor a. The sequence of phases reported is rich. For triangular lattices mostly, and although we were only 
interested in the lowest energy configuration, we have uncovered the relevance of the flexibility of the n-mers in a 
given trap, whereas previous approaches envisioned those composite objects as a rigid entity. In most cases, the 
qualitative predictions derived within the rigid picture appear qualitatively correct though. We have also addressed 
the adequacy of truncating interactions to colloids located in nearest neighbour traps, and found that in some cases, 
the quasi degeneracy between states of different orientational orders makes second nearest neighbours and possibly 
more remote shells relevant. 

We emphasise that the pinning potential we considered is of isotropic form, and therefore does not favour any 
orientation of the trapped n-mers. On the other hand, it was assumed in Ref. [Tl| that the pinning potential itself 
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leads to a discrete set of possible orientations, which were supposed to match the lattice symmetry. In particular, 
such an assumption would only be compatible with the ferromagnetic Po,o an d antiferromagnetic -Po.tt/ 2 phases for 
dimers on the square lattice. Our investigation shows that several other tilted phases exist, with an angular selection 
due to Coulombic repulsive interaction alone. Instead of a discrete set of predefined angles we find that the angular 
coordinates change continuously with parameters. We have explored in detail the limits of validity of both assumptions 
and have concluded that working with discrete angles can only be justified if the light confinement is very strong or 
when colloids repel weakly. We have shown however that the angular confinement is very weak in this regime. 
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